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Abstract. On the micro- and nanoscale, classical hydrodynamic boundary conditions such as the no-slip 
condition no longer apply. Instead, the flow profiles exhibit "slip" at the surface, which is characterized by 
a finite slip length (partial slip). We present a new, systematic way of implementing partial-slip boundary 
conditions with arbitrary slip length in coarse-grained computer simulations. The main idea is to represent 
the complex microscopic interface structure by a spatially varying effective viscous force. An analytical 
equation for the resulting slip length can be derived for planar and for curved surfaces. The comparison 
with computer simulations of a DPD (dissipative particle dynamics) fiuid shows that this expression is 
valid from full-slip to no-slip. 
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1 Introduction 

In the last years much effort has been spent on the de- 
velopment of new microfluidic devices. Recent success in 
reaching the nanometer scale has stimulated rising inter- 
est in fluid mechanics on submicrometer scales. Coarse- 
grained computer simulations are powerful theoretical tools 
for investigating the dynamical properties of such con- 
fined flows. Simulation studies rely heavily on the ade- 
quate choice and implementation of the boundary con- 
ditions at the surfaces. Traditionally, solid surfaces are 
taken to be "no-slip" boundaries, i.e., the fluid particles 
close to the surfaces are assumed to be at rest relative to 
the surface. While this can be rationalized to some extent 
for macroscopic (rough) surfaces [1] , it has no microscopic 
justiflcation. In fact, already J. C. Maxwell argued on sim- 
ple kinetic grounds that there always has to be some slip 
at surfaces, i.e., the tangential velocity does not entirely 
vanish at the boundaries [2 3 . In recent years, experiments 
have indicated that the no-slip boundary condition is in- 
deed usually not valid on the micrometer scale [4|5j . In- 
stead, it has to be replaced by the "partial-slip" boundary 
condition 

Sb 9nW|||rB = W|t|rB, (1) 

where denotes the tangential component of the veloc- 
ity and i9ni'|| its spatial derivative normal to the surface, 
both evaluated at the position of the so-called "hy- 
drodynamic boundary" . This boundary condition is char- 
acterized by two effective parameters, namely (i) the slip 
length Sb and (ii) the hydrodynamic boundary r^. We 
note that the latter does not necessarily coincide with the 



physical boundary. A boundary condition of the form ([T]) 
can be formulated very generally, even in situations where 
the "physical boundary" is not well-defined, e.g., on sur- 
faces covered with a coating or a wetting layer. Eq. ([1]) has 
been shown to be the most general hydrodynamic bound- 
ary condition compatible with conservation laws and irre- 
versible thermodynamics p]. 

The question arises how to implement such partial- 
slip boundary conditions in computer simulations. De- 
spite the experimental observation, most coarse-grained 
simulation models have aspired to realize no-slip bound- 
ary conditions. This task already turned out to be much 
more challenging than one might expect. For the dissi- 
pative particle dynamics (DPD) simulation method [6l7j . 
several papers have been devoted specifically to the prob- 
lem of reaching no-slip at the boundaries. Solid boundaries 
have been modeled, e.g., by frozen fluid regions [5t5] . or by 
planar walls containing one or several layers of embed- 
ded particles [9ll0lllj . It is worth noting that the presence 
of frozen embedded particles alone does not warrant no- 
slip ~ they have to be supplemented by a bounce-back 
reflection law at the walls. However, these may generate 
artifacts in the temperature profiles [9]. Approaches to im- 
plementing arbitrary partial-slip boundary conditions in 
Lattice-Boltzmann have recently been proposed by Benzi 
et al. |12I13| . However, to our best knowledge, no such 
method has been formulated yet for coarse-grained parti- 
cle simulations. 

In this paper, we propose an alternative way of im- 
plementing hydrodynamic boundary conditions in general 
coarse-grained simulations (DPD and other). Compared 
to the above mentioned methods, it has several advan- 
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tages. (i) It does not use embedded or frozen particles; 
hence, it does not introduce an artificial lateral structure 
in the walls, and it is cheaper from a computational point 
of view, (ii) It does not require the unphysical bounce- 
back reflections and thus avoids their possible artifacts, 
(iii) The basic idea is very general; it can be implemented 
for a wide range of simulation methods in a very straight- 
forward way. (iv) Most importantly: The method gener- 
ates full-slip, partial-slip as well as no-slip boundary condi- 
tions, and it provides a model parameter a with which the 
slip length can be tuned continuously and systematically 
from full-slip to no-slip. Moreover, an analytic equation 
can be derived, which allows to calculate the slip length 
to a very good approximation as a function of a. 

2 The basic idea 

The hydrodynamic boundary conditions at the surface 
result from interactions between the fluid particles and 
the walls. Depending on the microscopic structure of the 
wall/fluid interface, these can be quite complex. In our 
approach, we replace the boundaries by hard planar sur- 
faces, and the unknown atomistic forces by an effective 
coarse-grained friction force between the fluid particles 
and the walls. This leads to a dissipation of the kinetic 
energy and therefore to a decelerated velocity of the fluid 
close to the boundaries. The resulting slip length depends 
on the strength of the friction force. At friction force zero, 
one has a full-slip surface, i.e., slip length infinity. With 
increasing friction strength, the slip length decreases. 

The idea can be applied very generally. In the follow- 
ing, we will focus on particle-based off-lattice simulations. 
The friction force can then be implemented by introduc- 
ing spatially varying Langevin forces acting on particles 

Ff = Ff + Ff (2) 
The dissipative contribution 

Ff = -7L UJL{zi/Zc) {Vi - Vumll) (3) 

couples to the relative velocity (v^ — Vwaii ) of the particle 
with respect to the wall, with a locally varying viscosity 
jLLOLizi/ Zc) that depends on the wall-particle distance Zi, 
with a cutoff distance Zc- The weighting function ll!l{t) is 
positive for T < 1 and zero for t > 1. Otherwise, it can 
be chosen freely. The prefactor can be used to tune 
the strength of the friction force and hence the value of 
the slip length. To preserve the global temperature T and 
to ensure the correct equilibrium distribution, a random 
force obeying the fluctuation-dissipation relation has to 
be added, 

F«„ = V/27L kBTLJL{z^/Zc) X^,a (4) 

with a — x,y, z, where fc^ is the Boltzmann constant and 
Xi,a a Gaussian distributed random variable with mean 
zero and unit variance: {xi,a) = 0, {Xi,aXj,0) = ^ij^ap- 
Eq. ^ can be used to model interactions with immobile 
walls (such as channel boundaries) as well as interactions 



with surfaces of mobile and/or rotating objects (such as 
colloids). In the latter case, the force ^ must be balanced 
by a countcrforce — Ff and a countertorque — (ri — R) x Fp 
acting on the object (R being its center of mass). 

An analytical expression for the slip length shall be 
derived in the next section. In short, we find that the 
result can be written as a function of the dimensionless 
parameter 

zl -iLPh, (5) 

where p is the fluid density and r/ its shear viscosity. To 
leading order in a, the slip length is given by 

— = (6) 

Zc a dr w(t) 

The sign of the next-to-leading correction is usually neg- 
ative. For example, a steplike weight proflle lo[t) = 1 for 
T < 1, a;(r) = otherwise, gives 

-55 1 2 1 
-=a-3~45" + ---- 

(cf. eq. ^9]) and a linear function ujir) — 1 — r 



Zc a 15 1800 ' ' 

(cf. eq. (HOI))- Hence the slip length becomes zero for an 
appropriate choice of a, and it can be tuned to any value 
up to infinity (corresponding to a = 0) . 



3 Theory 

The theory is based on one crucial assumption: The Navier- 
Stokes equations are taken to be valid on the length scale 
Zc of the cutoff. This assumption may seem bold, given 
that Zc will typically be chosen of the order of one parti- 
cle diameter. However, computer simulations .14,15j have 
shown that the Navier-Stokes equations are valid down to 
surprisingly small length scales. This is also supported by 
our own simulations (see below). 

For simplicity, we begin with discussing planar sur- 
faces. The coordinate system is chosen such that the sur- 
face is at rest and located at z = 0. The stationary Stokes 
equation for our system reads 

riv"(z) ^ pjL UJl{z/Zc) v{z) - p fext (9) 

with the fluid density p and the shear viscosity rj. Here, 
we assume that the fluid viscosity does not change in the 
vicinity of the walls. The external force fext incorporates, 
e.g., the effect of a pressure gradient in the direction of 
flow. The physical wall in our system is smooth, and we 
have no explicit solid/fluid wall friction. Thus the shear 
stress in the fluid must vanish, and the velocity field v{z) 
satisfies the boundary condition 

v'{z^O)^0. (10) 
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We recall that the viscous force lol{z / Zc) v{z) is only 
active within a layer of finite thickness Zc- The total fric- 
tion force per surface area generated in this layer is given 

by 

FIA^ p-iL r Azujl{z/z,)v{z). (11) 

Our goal is to replace the layer by an effective boundary 
between a solid and an unperturbed fluid, i.e., a hypothet- 
ical fluid not subjected to the additional viscous forces in 
the layer. The friction force per surface area, F/A, is then 
equal to the frictional stress on the solid, and to the shear 
stress on the unperturbed fluid at the position z = zb oi 
the boundary. This yields the effective boundary condi- 
tions 

^F/A^CBV^'>\zB)^nv^''^'{zB), (12) 

where is the unperturbed velocity profile, and C,b the 
fluid-solid friction coefficient. The comparison with eq. ([1]) 
shows that the slip length can be identified with 5b — 

v/Cb- 

We first show that the hydrodynamic boundary must 
be identical with the physical boundary, 

ZB = 0. (13) 

To this end, we first insert eq. ([9]) in eq. (fTT|) and perform 
a partial integration, taking advantage of the boundary 
condition (jlOp. to obtain 

- F/A = 7^v'{Zc)+ pfe^tZc- (14) 

According to eq. (fT2l) . this must be equal to rjv'''^^' {zb)- 
The unperturbed velocity profile v^^^ solves the Stokes 
equation at 7^ = with the boundary conditions 

«(0)(z,) = «(ze); «(0)'(z,)=t;'(z,). (15) 

In the planar case (HI), one obtains 

(z) = v{z,) + v'{z,){z - z,) - ^fe.t{z - z,f (16) 
and hence 

riv^°'''{zB) = riv{zc) + fextp{zc - zb)- (17) 

Comparing eqs. (fT7|) and (fT4|) gives eq. (fT3|) . 

Next we calculate the slip length. This is most con- 
veniently done for the case f^xt — 0. We first integrate 
eq. dHl) for given u(0) =: and w'(0) = (eq. dTO])) to 
get v[zc) and v' [zc). Via eq. (fTB|) . we then determine the 
values v^^^[zb) and v^^^'{zb) of the unperturbed profile at 
the hydrodynamic boundary. This finally allows to calcu- 
late the slip length from 

<5s=^;(°)(zB)/t;(")'(zs) (18) 

(cf. eq. ((T2)) V We note that all profiles scale linearly with 
vq, hence the final expression for the slip length does not 
depend on any more. More generally, the result depends 
only on the dimensionless quantity a defined in eq. ((S]), 




Fig. 1. Two types of cylindrically curved geometries. In the 
Couette geometry (left), the flow profile bends along the curved 
surface. In the parallel geometry, the fluid flows in the perpen- 
dicular direction. 



in the sense that the slip length in units of Zc, i.e., the 
quantity 5b/ Zc, can be written as a function of a only. 
For example, if the weight function wl(t) is a simple step 
function, lol{t) = 1 for r < 1, the flow profile v{z) at 
z < Zc is a superposition of exponentials exp(±Y^z), and 
one obtains 



Zc v^tanh(Q;) 

If drops linearly to zero, u}l{t) — 1 — t, the profile 

can be written as (see appendix A) 




where F is the Gamma- Function and / the modified Bessel 
function of the first kind. 

The slip length can become zero (corresponding to no- 
slip) or even negative. The no-slip boundary condition is 
obtained at a = 1.433 for a steplike weight function, and 
at a = 3.973 for a linear weight function. Negative slip 
lengths are encountered at even larger a. In this case, the 
hypothetical unperturbed profile v^°^ changes sign close to 
the boundary. We note that the true velocity profile, v{z), 
never changes sign; hence our negative slip lengths do not 
correspond to unphysical situations. 

In the regime of positive slip lengths, a is small and can 
be used as an expansion parameter. Expanding eqs. pop 
and (|20)) in powers of a gives ([7]) and ^ . More generally, 
one can derive a useful expression for arbitrary weight 
functions. We first note that the true velocity profile and 
the unperturbed profile are identical at the order a", i.e., 
v{z) — w^°)(z) -I- 0{a). Furthermore, the derivative of the 
velocity profiles is of order a, by virtue of eqs. (fTTI) . (fT^ . 
and P?)l . hence one even has v{z) = i;(°'(zb) -|- 0{a). 
Applying once more eqs. pT|) . (I12p . and we obtain 

eq. ([6]). This equation allows to estimate the slip length 
reasonably accurately for arbitrary choices of the weight 
function w(t). 

The extension of the theory to curved geometries is 
straightforward. Here we discuss cylindrical curvature (i.e., 
edges). One has to distinguish between curvature in the 
direction of flow (the "Couette" case) and curvature in 
the transverse direction (the "parallel" case), see Fig. [T] 
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The curvature is characterized by the radius i? of a tan- 
gent sphere, which we take by convention to be positive at 
the surface of convex objects, and negative at the surface 
of concave objects. Details of the calculation are given in 
the appendix B. As in the planar case, the hydrodynamic 
boundary is found to be identical with the physical bound- 
ary. To leading order of a, the slip length is given by 



dTa;(r) (l + r-^)2-f 0(a"). (21) 
K 



in the Couette case, and 



Zc_ 

Sb 



drw(r) (1 



r|)+0(a°). 



(22) 



in the parallel case. In the Couette case, the inverse slip 
length has a contribution of purely geometric origin. 
The existence of this term has been pointed out by Einzel 
et al. [1]. The remaining term can be identified with a 
"microscopic" inverse slip length. Both in the Couette and 
the parallel case, this microscopic slip length is reduced at 
the surface of convex objects, and enhanced at the surface 
of concave objects. 



4 Comparison with computer simulations 

To test the method and the theory, we have carried out 
dissipative particle dynamics [617] (DPD) simulations of 
fluids confined in a parallel slit. The fluid particles in- 
teract with purely dissipative DPD forces of range a 7 . 
They have no conservative interactions with each other, 
their static equilibrium structure is thus that of an ideal 
gas. They are confined into the slit by two parallel walls, 
with which they interact via a Weeks-Chandler-Andersen 
potential |16| with characteristic length a; alternatively, 
one can use hard reflecting walls. In the vicinity of the 
walls, up to the cutoff distance Zc, the particles experi- 
ence a wall friction force of the form ([2]) with a linear 
weight function, ujl(t) = 1 — t. The choice of the weight 
function is in fact completely arbitrary. It was motivated 
by the fact that the DPD weight factors for intermolec- 
ular dissipative interactions are also usually chosen lin- 
ear for reasons of computational efficiency. The natural 
units of the simulation are the length unit tr, the ther- 
mal energy ksT =: e, and the mass of the particles m. 
In these units, the cutoff at the walls was chosen Zc — cr, 
and the time step 5t = 0.01 a^m/e. Periodic boundary 
conditions were applied in the slit plane. All simulations 
were carried out with extensions of the simulation package 
ESPResSo [17118119]. 

We have simulated both planar Couette and planar 
Poiseuille flow [20]. In the flrst case, one of the walls moves 
at constant speed relative to the other. In the second 
case, the effect of the pressure gradient is mimicked by 
a constant force acting on all particles in the slit. Fig. [2] 
shows two examples of flow profiles. Sufficiently far from 
the walls, the classical parabolic Poiseuille flow and linear 
Couette flow profiles are nicely recovered. 
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Fig. 2. Normalised plane Poiseuille flow (top) and plane Cou- 
ette flow (bottom) for DPD fluids {^dpd ~ S^/me/a) with 
number density p = 3.75a^^. The parameters of the surface 
friction force are = 1.0^/me/a, Zc — a. The dashed lines 
show fits to the theoretical profiles. 



The combined results from the Poiseuille and Couette 
fits allows to calculate the slip length and the position 
of the hydrodynamic boundary independently: We deter- 
mine the distance P between the two points where the ex- 
trapolated Poiseuille parabola vanishes, and the distance 
C between the two points where the extrapolated Couette 
line reaches the velocities of the two walls. The slip length 
is then given by 



4 = (C2 - p2)/4^ 



(23) 



and the hydrodynamic boundary is located at the distance 



ZB = &B-{C- L)/2. 



(24) 



from the physical boundary, where L is the distance be- 
tween the two physical walls. We note that both C and 
P enter the expression for the slip length, hence both the 
Poiseuille and the Couette profile have to be simulated. 
This is because we do not make any assumptions regarding 
the hydrodynamic boundary. If its position is unknown, 
two types of profiles have to be measured in order to de- 
termine the two parameters zb and 5b- This fact is not 
always appreciated in the literature. 

We have determined the hydrodynamic boundary and 
the slip length for the range of parameters 7l — (0.1 - 
5) ^me/a (the friction at the wall), p = (3.75-12.5) 
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Fig. 3. Position of the hydrodynamic boundaries for varying 
yoPD and fixed number density p — 3.75a~^ (top) and for 
varying rho and fixed jdpd ~ 2y/me/a (bottom). The result 
does not depend on the value of the surface friction 7_l. The 
data shown here are averaged over all 7^. 



(the number density of fluid particles), and 'Jdpd = (1- 
10)-y/me/cr (the friction coefficient of the DPD forces). The 
results are summarized in the figures [3] and 21 The the- 
ory predicts that the hydrodynamic boundary is identical 
with the physical boundary, eq. (fTS]) . This was found to be 
correct for most systems. Significant deviations were only 
encountered at small DPD-friction -foPD (fig- 131) • Fig. 2] 
shows the results for the slip length as a function of the di- 
mensionless parameter a = jLPZc/v- The shear viscosity 
r] of the fluids was determined from the Poiseuille fitsf^. 
The numerical results for the slip length compare very 
favorably with the theoretical prediction of eq. ([8]), ex- 
cept in the regime where the hydrodynamic boundary also 
deviates from its theoretical position, i.e., at small DPD- 
friction "fDPD- Even in this regime, the method can still be 
used to produce partial slip. A wide range of slip lengths 
has been realized. The inset of Fig. 2] illustrates that the 
no-slip regime can also be reached very naturally. 
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Fig. 4. Slip length 5b in units of Zc vs. a for varying val- 
ues of the parameter triplet (p,^DPD,'yL) (in units of 
or ^yme/a, respectively). Black: series with p fixed: (3.75 , 2- 
10, 0.1-1). Red: series with jdpd fixed: (3.75-12.5, 2, 0.1-1). 
Blue: selected triplet values: (6.35,5,1), (5, 5,1), (11.25,2,1.1), 
(11.25,2,1.2), (3.75,10,2.5). Green star symbols: series corre- 
sponding to the point in fig. O where the position of the 
hydrodynamic boundary deviates from theoretical expecta- 
tion: (3.75,1,0.1-1). Dashed black line: Theoretical prediction 
of eq. (|20|l . The inset shows a blowup of the same data. 

5 Summary 

We have presented a method that allows to implement 
arbitrary hydrodynamic boundary conditions in coarse- 
grained simulations. We have illustrated and tested this 
method with DPD simulations; however it can be applied 
much more generally. The main idea is to introduce an 
effective surface friction force, complemented by an ap- 
propriate thermostat. This idea can also be implemented, 
e.g., in stochastic rotation dynamics simulations or in Lat- 
tice Boltzmann simulations. The strength of the friction 
force effectively determines the slip length. We have shown 
that our method allows one to tune the slip length system- 
atically from full-slip to no-slip. Furthermore, we have de- 
rived an analytical equation for the slip length as a func- 
tion of the strength of the friction force and the shear 
viscosity of the fluid. 

An alternative force-driven method to control the slip 
length on surfaces has recently been proposed for the Lat- 
tice Boltzmann simulation method by Benzi et al[12j. In 
this work, conservative forces are introduced that deplete 
the fluid density at the surface, which leads to partial slip 
- a physical effect. In contrast, the dissipative forces in our 
method mimick directly the effective friction between the 
surface and the bulk fluid. While the method does thus 
not contribute to a microscopic understanding of the slip 
at surfaces, it does offer a way to model the essential wall 
properties, from the hydrodynamics point of view, and 
to absorb the complex wall structure into a well-defined 
coarse-grained parameter set [13]. The method can be im- 
plemented easily, it generates the partial-slip boundary 
conditions for microfluidic simulations that are observed 
in experiments. Curvature can be introduced in a straight- 
forward manner, the effect of curvature on the slip length 
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in the model can be calculated analytically, and corrected 
by an easy change of local friction if necessary. Hence the 
method is suitable to study complex geometries and/or 
objects, such as structured channels, channels with ob- 
jects, or rotating nanoparticles in flow. 

The only obvious restriction is that the characteris- 
tic length scale of the structures has to be larger than the 
length scale on which the friction force is applied. In fu- 
ture work, it will be interesting to determine the minimal 
lateral length scale that the method can handle, t.g.^ by 
studying surfaces with mixed boundary conditions. For ex- 
ample, analytical results are available for the macroscopic 
effective slip length on surfaces covered by stripes with al- 
ternating full-slip and no-slip boundary conditions |22I23| . 
and the Stokes equation has also been solved numerically 
for more general mixed surfaces [24125] . This will be a good 
starting point for a systematic study. Furthermore, we 
plan to apply our method to the simulation of electrolyte 
flow to assess the influence of partial-slip boundary con- 
ditions on electroosmotic flow 26 . 

We thank Burkhard Diinweg, Christian Holm and Ulf 
D. Schiller for nice and fruitful discussions. This work was 
funded by the Volkswagen Stiftung. The visit of MPA to 
Bielefeld was supported by the Alexander von Humboldt 
foundation. 



6 Appendices 

A: Slip length for linear weight function 

The program for calculating the slip length (Eq. (fT8|l can 
easily be carried out numerically for arbitrary weight func- 
tions w(t). For some simple choices of the weight function, 
analytical results can also be obtained. To illustrate the 
approach, we sketch the calculation leading to Eq. ([^0]) . 
the expression for the slip length for a linear weight func- 
tion lo(t) — \~t. 

In the absence of an external force (/ext = 0), the 
general solution of Eq. ([9]) for a linear weight function 
reads 

v{z) = A Ai(c(l - z/zc)) + B Bi(c(l - zjz^)) (A.2) 

with c — — (— a)^/^, where Ai and Bi are the Airy func- 
tions and a has been defined in Eq. ([5]). The boundary 
conditions w(0) = wq: ^^'(0) = determine the coefficients 

A = woAi'(c)/(Ai'(c)Bi(c) - Ai(c)Bi'(c)) 
B = woBi'(c)/(Ai'(c)Bi(c) - Ai(c)Bi'(c)). 

The unperturbed profile is linear, u(°^(z) — w(zc)+'y'(zc)(2- 
2c) 7 and the hydrodynamic boundary is located at = 
according to Eq. , hence the equation for the slip 
length can be written as 



Bi(0) = V3Ai(0), and Bi'(O) = -V3Ai'(0), where F is the 
Euler gamma function, one obtains 

'^B^ 1 Ai-(c)73-Bi'(c) 

3V3c Ai'(c)V3 + Bi'(c)' ^ ■ ' 

To get rid of the complex argument c, we use the series 
expansion of the Airy functions [2 7[. The functions of in- 
terest for us can be written as Ai'(z)\/3 = 1\ — I2 and 
Bi'(z) = /i +/2 with 



(A.3) 



Inserting this in Eqs. (jA.2p with (|A.3p and (|A.3p and using 
the identities Ai(0) = l/(32/3r(|)), Ai'(O) ^ -l/(3i/3r(i 



331/6 ^ rik + hli) k\ ' 9 

k 



_3 



r{k + l/3) kl ' 9 



Comparing this with the series representation of the mod- 
ified Bessel function of the first kind. 



2' ^ r{k + iy + l)k\ '2 



(A.5) 



one easily verifies the identities 



2V^ ^2/3 2^ 

/i=-^/-V3(^) and /2 = ^%3(^). 

Inserting this into Eq. (|A.4|) gives Eq. P0|). 



B: Curved boundaries 

The strategy for calculating the slip length and the hydro- 
dynamic boundary position at curved surfaces is perfectly 
analogous to that sketched in the main text for the planar 
geometry. We begin with noting that the expression for 
the total friction force per surface area, eq. pT|) . acquires 
an additional geometric factor at curved surfaces. 



F/A = pjL 



ArLOL{{r-R)/z,) v{r)—, (A.6) 



which accounts for the fact that the viscous layer is stretched 
or compressed in the presence of curvature. Here, R is the 
radius of curvature, r the distance from the center of the 
tangent sphere, and tb the position of the hydrodynamic 
boundary. Furthermore, the Stokes equation, ([5]), and the 
boundary conditions, (fTO|) and (fT2|) have to be replaced by 
the appropriate cylindrical versions. 

We first discuss the Couette case. The Stokes equation 
for this geometry reads 

■qdr-drrv{r) = p7^a;((r - R)/zc) v{r) - p fext, (A. 7) 
r 

The unperturbed solution of this equation for 7l = and 
fixed v{rc) = Vc,v'(rc) = v'^ is given by 



V 



'"\r) = i[f (-c + .>c + /) + ^(.e-«^r.-{)-{(f )2] 



(A., 
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with / = pfexti'c/'n- I'^ the absence of an external force, 
f — 0, one recovers the well-known profile of cylindrical 
Couette flow. The appropriate expressions for the fluid 
shear stress also differ from those in the planar case, and 
Eqs. Uni) and ([12]) have to be replaced by[T] 



{drV - v/r)\r=B. = 0, 



(A.9) 



- F/A = CBV^°\r)\r^rn = V{drV^°^ " V^°^ / r)\r^rs ■ 

(ATO) 

Comparing cq. (jA.lOp with cq. ([T]), one immediately ob- 
tains 6g = Sq + with 6q = Cb/V: where is a 
geometric contribution, and 5^^ incorporates the specific 
effect of the surface structure. Comparing with eq. (IA.6|) . 
one gets 



- = a dr UJL[T) — 
Jo V 



v{R + TZc) R + TZc 



(A.ll) 



For fixed v{R + Zc) = Wc, v'{R + Zc) = v'^, the actual pro- 
file v[r) is identical with the unperturbed profile v^^\r) of 
eq. (|A.8|) at leading order of a. Inserting this in eq. (jA.lip 
and using tb = R, we finally recover eq. ((2T|) . To establish 
the identity tb — R, we require again that the total fric- 
tion force as given by eq. (|A.6[) is identical with the hypo- 
thetical shear stress on the unperturbed fluid, eq. (jA.lOp . 
Inserting (|A.7j) in (jA.6p . carrying out a few partial inte- 
grations, using (jA.9[) . and finally replacing once more v{r) 
by -y'^'^^(r) at order a, we obtain 



rc 



Comparing this with eq. (jA.lOp . after inserting (jA.Sp . we 
find that tb must be equal to the surface position R at 
the leading order of a. 

The parallel case is even simpler. The appropriate Stokes 
equation reads 

rj-drrdrv{r) = p-fLLo{{r - R)/zc) v{r) - p fext- (A. 13) 



In the unperturbed case = 0, this is solved by 

«(«)(r) =v, + + |)ln(-) + |(1 - (-)') (A.14) 
2 Tc 4 rc 

with Vc := v{rc), v'^ := v'{rc), and / = pfextrl/rjas before. 
The expressions for the fluid shear stress are similar to 
those in the planar case. 



drV\r=R = 0, 



(A.15) 



References 

1. Einzel D., Panzer P., and Liu M., Phys. Rev. Lett. 64, 
(1990) 2269. 

2. MaxweU J. C, Philos. Trans. Soc. London 170, (1867) 
231. 

3. Kennard E. H., Kinetic theory of gases (McGraw-Hill, 
New York, 1938). 

4. Pit R., Hervet H., and Leger L., Phys. Rev. Lett. 85, 
(2000) 980. 

5. Neto C, Evans D. R., Bonaccurso E., Butt H.-J., and Craig 
V. S. J., Rep. Prog. Phys. 68, (2005) 2859. 

6. Hoogerbrugge P. J. and Koelman J. M., Europhys. Lett. 
19, (1992) 155. 

7. Espafiol P. and Warren P., Europhys. Lett. 30, (1995) 
191. 

8. Boek E., Coveney P., Lekkerkerker H., and van der Schoot 
P., Phys. Rev.E 55, (1997) 3124. 

9. Revenga M., Zuniga L, and Espafiol P., 
Comp. Phys. Comm. 121-122, (1999) 309. 

10. Duong-Hong D., Phan-Thien N., and Fan X. J., Compu- 
tational mechanics 35, (2004) 24. 

11. Pivkin L V. and Karniadakis G. E., J. Comput. Phys. 
207, (2005) 114. 

12. Benzi R., Biferale L., Sbragaglia M., Succi S., and Toschi 
F., Europhys. Lett. 74, (2006) 651. 

13. Benzi R., Biferale L., Sbragaglia M., Succi S., and Toschi 
F., J. Fluid Mech. 548, (2006) 257. 

14. Alder B. J. and Wainwright T. E., Phys. Rev. A 1, (1970) 
18. 

15. Diinweg B. and Kremer K., J. Chem. Phys. 99, (1993) 
6983. 

16. Weeks J. D., Chandler D., and Andersen H. C, J. Chem. 
Phys. 54, (1971) 5237.^ 

17. ESPResSo-Homepage, |http://www.espresso.mpg.de| 

18. Arnold A., Mann B. A., Limbach H.-J. and Holm C, 
GWDG-Bericht: Forschung und wissenschaftliches Rech- 
nen 2003 Vol. 63, Kremer K. Macho V. Eds. (Gesellschaft 
fiir wissenschaftliche Datenverarbeitung mbh, Gottingen, 
Germany, 2004), pp. 43 ff. 

19. Limbach H.-J., Arnold A., Mann B. A., and Holm C, 
Comp. Phys. Comm. 174, (2006) 704. 

20. Lu P.-C, Introduction to the mechanics of viscous fluids 
(Holt, Rinehart and Winston Inc., New York 1973). 

21. Backer J. A., Lowe C. P., Hoefsloot H. C. J., and ledema 
P. D., J. Chem. Phys. 122, (2005) 154503. 

22. Philip J. R., J. Appl. Math. Phys. (ZAMP) 23, (1972) 
353. 

23. Lauga E. and Stone H.A., J. Fluid Mech. 489, (2003) 55. 

24. Cottin-Bizonne C, Barentin C, Charlaix E., Bocquet L., 
and Barrat J.-L., Eur. Phys. J E 15, (2004) 427. 

25. Priezjev N.V., Darhuber A.A., and Troian S.M., Phys. 
Rev. E. 71, (2005) 041608. 

26. Smiatek J. et al. , Poster presented at the International 
Soft Matter Conference, Aachen, October 2007. 

27. See, e.g., functions.wolfrain.com. 



- F/A = <:BV^'>Hr)\r=rn = VidrV^^^ " V^''> / r)\r=rn ■ 

(A.16) 

Proceeding as before, we find that = i? is fulfilled ex- 
actly, and we recover eq. (|22l) . 



